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Abstract 

We study regular and black hole solutions to the coupled classical Einstein- Yang-Mills- 
Higgs system. It has long been thought that black hole solutions in the spontaneously broken 
phase of such a theory could have no nontrivial field structure outside of the horizon. We 
first show that the standard black hole no-hair theorem underlying this belief, although true 
in the abelian setting, does not necessarily extend to the non-abelian case. This indicates 
the possibility of solutions with non-trivial gauge and Higgs configurations decaying expo- 
nentially outside the horizon. We then find such solutions by numerical integration of the 
classical equations for the case of SU(2) coupled to a Higgs doublet (the standard model 
less hypercharge). As a prelude to this work we also study regular and black hole solutions 
to Einstein-Non-Abelian-Proca theory and as a postscript we briefly discuss the important 
issue of stability. 



I. Introduction 



Some time ago it was proven that black hole solutions in Einstein gravity coupled to the 
abelian Proca model or to a spontaneously broken abelian gauge theory cannot have any 
nontrivial field structure outside of the horizon [1-2]. In conjunction with similar arguments 
applied to a variety of field theories coupled to gravity [1-4] , a widely held belief has been 
that the only distinguishing features of a black hole exterior to the horizon are charges 
carried by massless gauge fields [6]. These are mass, angular momentum and "electric" or 
"magnetic" charges, where the latter refer to charges carried by some possibly non-abelian 
gauge group. This belief is often referred to as the black hole no-hair theorem (or more 
precisely, no-hair conjecture). 

In elementary particle physics we are familiar with the fact that non-abelian symmetries 
often play a crucial role. Based upon the abelian case and more general prejudice, it has 
generally been accepted that the no-hair statements ensure that no evidence of the broken 
classical charges can be found outside a black hole horizon. In particular, it has been 
thought that one cannot find stationary black hole solutions which are nonsingular at the 
event horizon and have a nonvanishing massive vector field that decays exponentially outside 
of the horizon (for example, see [1-2], [5] and [7-8]). It is the intent of the present paper to 
quantitatively investigate the veracity of this belief. That is, we study black hole solutions 
in Einstein gravity coupled to the non-abelian Proca model and to a spontaneously broken 
gauge theory (with symmetry group SU(2) in each case) in order to determine whether 
the no-hair statements are affected by the non-abelian structure. Our results indicate that, 
in fact, the non-abelian nature of the theory has a crucial impact giving rise to classical 
solutions which violate the no-hair conjecture. 

In particular, we numerically study spherically symmetric solutions to Einstein gravity 
coupled to SU(2) gauge theory and a Higgs doublet, as in the standard model (without 
hypercharge) . We take the standard model form for the symmetry breaking Higgs potential 
with only the lower real component of the Higgs field having a nonzero vacuum expectation 
value. Integrating via a two parameter shoot (associated with the initial values of the gauge 
and Higgs fields), we find black hole solutions in which there are nontrivial gauge and Higgs 
fields outside of the horizon (with no global charges associated with the gauge or Higgs fields) 
that ultimately exponentially decay to their vacuum values. In fact, we find two families 
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of such solutions with each member distinguished by the number k of zero crossings of the 
single free function in the gauge connection. The k — 1 member of one of these two families 
of solutions has a limit (the dimensionless Higgs vacuum expectation value going to zero) 
in which gravity becomes arbitrarily weak outside of the horizon and the solution resembles 
the known flat space sphaleron solution [9]. So, at least this particular solution might be 
characterized as a being a black hole with sphaleron hair. Our study indicates, however, 
that the other solutions do not have a weak gravity limit and hence do not have known flat 
space counterparts. 

We note that ours is not the first challenge to the no-hair conjecture, but it is useful 
and important to emphasize the viewpoint espoused in [5]. As these authors point out, 
one should distinguish primary hair from secondary hair. The latter refers to black hole 
structures which exist solely as the result of (well known) primary hair such as gauge charges 
and hence are not fundamentally new characteristics. The existence of secondary hair, 
therefore, is not really in conflict with the no-hair conjecture. The charged dilaton black 
hole solutions [10-15] in which there is nontrivial scalar field configuration outside of the 
horizon provides an example of this distinction. The dilaton configuration is nontrivial 
because the electric charge (primary hair) acts as a source and hence yields only secondary 
hair. Another interesting solution which (among other things) provides a challenge to the 
no-hair theorems are the "black holes inside magnetic monopoles" found in [16, 17]. Again, 
however, the hair associated with the nontrivial field configuration outside of the horizon 
is of the secondary sort [5]. Two other challenges to the no- hair conjectures, which are 
closer to having primary hair, are the works of [18] (motivated by the paper of Bartnik and 
McKinnon [19]) and [20]. The work of [18] finds numerical black hole solutions of Einstein 
gravity coupled to SU(2) gauge theory in which the fields decay sufficiently quickly so as to 
have vanishing global charges. The hair on these solutions is still carried by massless gauge 
fields - the novel feature is that the fields leave no imprint at infinity. Unfortunately, these 
solutions have been shown to be unstable [21]. In [20] some interesting numerical work also 
has shown that black holes can exist with nontrivial Skyrme fields in the Einstein-Skyrme 
model. 

Our work on this subject was partly motivated by the publication of [19] in which 
the authors found smooth solutions to the coupled Einstein-S'[/(2) gauge theory classical 
equations [23]. These smooth solutions were surprising in that no-go theorems for classical 
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glueball solutions to non-abelian gauge theory (without coupling to gravity) were established 
in [24] while no-go theorems for smooth solutions including gravity were proven for 2 + 1 
dimensions in [25]. Unfortunately, the 3 + 1 dimensional solutions of [19] were subsequently 
proven to be unstable. A natural question to ask is whether such smooth solutions in 
3 + 1 dimensions can be found in the physically relevant case of spontaneously broken non- 
abelian gauge theories and, if so, whether they might be stable. In addition to our black 
hole solutions, we do in fact find and present such smooth solutions but have not as yet 
exhaustively analyzed their stability properties. 

Stability is, in fact, an important unanswered question about both our black hole and 
regular solutions. However, we hasten to emphasize that the long held belief that such black 
hole solutions do not exist is based on the no-hair results which themselves have nothing to 
do with stability. Rather, the no-hair results were based on the careful examination of certain 
classical field theories coupled to gravity which simply have a different character from the 
non-abelian vector theories studied here. However, as far as physical relevance is concerned, 
stability is a crucial feature. General prejudice would certainly lean towards suspecting that 
our solutions are unstable for two main reasons. First, without a global charge leaving an 
imprint at infinity, stability certainly seems less likely. Second, as mentioned, there is a 
relation between our solutions and flat space sphalerons, the latter of which are unstable. 
On the other hand, some of our solutions seem to involve the trapping of gauge bosons in a 
region of space where, due to the Higgs configuration, they are less massive (relative to their 
asymptotic mass), and this, potentially, may bode well for stability. It has been noted that 
there might be a no-hair theorem if one only considers stable solutions [21]. We do not know, 
however, of any definitive arguments supporting this natural suggestion. We are presently 
studying the stability question for both the smooth and black hole solutions discussed here 
and will report on this elsewhere. 

The organization of this paper is as follows. In section II we review the no-hair the- 
orems and indicate why they do not apply in the situation under study. This opens up 
the possibility for the existence of black hole solutions with non-abelian gauge and Higgs 
hair. In section III we embark on the numerical construction of such solutions, beginning 
with the simpler Einstein-Non-Abelian-Proca system for the case of SU(2). This system 
is interesting in its own right as a first concrete illustration of how the classical solitons in 
abelian and non-abelian Proca theories differ substantially. It also is a useful prelude to 
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our discussion of the Einstein- Yang-Mills-Higgs system in section IV. In both sections III 
and IV we present spherically symmetric regular and black hole solutions. In section V we 
present some preliminary work on the stability and we offer our conclusions. 



II. Limitations on No-Hair Theorems 

Many of the no-hair proofs for classical field theories utilize an elegant and powerful 
method first developed by Bekenstein [1]. It is, in effect, a repackaging of the generally 
covariant field equations into a statement about the behavior of fields outside the event 
horizon. In this section, we briefly review Bekenstein's arguments and illustrate how non- 
abelian gauge fields in general, and the Einstein-Non-Abelian-Proca (ENAP) and Einstein- 
Yang-Mills-Higgs (EYMH) systems in particular, can possess static black hole solutions with 
nontrivial exterior structure. 

Consider the action for a set of arbitrary local fields $k in a gravitational background: 
S = J d 4 xC = J d 4 Xy/—gC. After multiplication by cftx&k and integration by parts, the 
covariant Euler-Lagrange equations become 
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The left-hand side of this equation can be expressed as a surface integral dS^, where 
dS/j, is an element of the hypersurface bounding the volume Q over which we integrate, and 
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(2.2) 



If we assume that our system admits black hole solutions and we choose our four-volume to 
be the black hole exterior, then dfl consists of the horizon, spatial infinity, and future and 
past timelike infinities. The behavior of at spatial infinity is already determined for all 
nontrivial and physically relevant fields by the field equations: — > 1/r 3 asymptotically for 
massless fields and 6 M vanishes exponentially for massive fields, so there is no contribution 
at spatial infinity to the left side of (2.1). Since dS^ is proportional to the normal n M of 
the hypersurface, which can be chosen to satisfy n; = as t — > oo, dS^ also vanishes 
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at timelike infinity when 6° = 0. This condition on b° is satisfied by static fields, and thus 
there are no non- horizon contributions to the surface integral in eq.(2.1). From the fact 
that the horizon is a null hypersurface [gijdS l dS J = 0) and that gij is positive semi-definite 
on the horizon, it can easily be shown that dS^ vanishes on the horizon when b^ is 
bounded there. Thus, for static fields with finite b^ on the horizon, eq.(2.1) implies that 
the integral over the entire black hole exterior of some nontrivial function of the fields and 
their derivatives must vanish. The strategy for establishing no hair becomes clear: if one 
can demonstrate that the integrand in eq.(2.1) is negative or positive definite, then the only 
finite energy solutions satisfying (2.1) are those for which the integrand vanishes. For most 
Lagrangians of physical interest, this implies that the fields must assume constant values 
over the entire black hole exterior. 

The Bekenstein approach has been used to establish no-hair theorems for Klein-Gordon 
theory, abelian Proca [1] and Higgs theories, and the Goldstone model [2], among others. 
Here we provide a sketch of the right-hand side of eqn.(2.1) for non- abelian gauge fields which 
demonstrates the possibility of the existence of Einstein- Yang-Mills (EYM) black holes. 

The Lagrangian for an SU (2) gauge theory may be written 



-EYM = 

4"7T 



where i is the isospin index. A straightforward calculation of the right side of eq.(2.1) for a 
metric with signature (—1, +1, +1, +1) gives 

__L J (-87r£ EYM + l - [ge ijk Af A^] F®^ , (2.4) 

(i) (i) (i) (j) (k) 

where g is the gauge coupling: F^ = d^At,' — d v Ay + ge^Aji 1 A\j . If we consider only 

(i) 1 2 

static fields, and for simplicity assume Ay =0, then |F| > 0. The second term in the 
integrand, however, is not necessarily positive: its sign will depend on the details of the 
solution in the black hole exterior Q, so it is clear that the nonlinear terms can lead to a 
possible way of avoiding this no-hair argument. 

To demonstrate this point more quantitatively, consider the expression in (2.4). Let 
Cp,l = d^A^ — d v A$ and = ge l - k A^ A^\ Then to avoid the Bekenstein argument 
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we require the non-positivity (somewhere) of S = (C + D) ■ (C + D) + D ■ (C + D) where 
C-C — C^ 1/ C^ fiV etc. With the gauge choice A*jp = 0, and the assumption that in this gauge 
all fields are time independent (no 'electric fields'), it follows that (aC + [3D) ■ (aC+/3D) > 
for any real a, (3. The relation S = (C + 2D) ■ (C + D) < requires that C ■ D < 0, and 
gives C ■ C + 2D ■ D < —3C ■ D. This may be rephrased as 

hc-C + D-D) - hc-C- D-D) < -C-D < -(C-C + D-D), (2.5) 

where the second inequality is just the positivity of (C + D) • (C + D). 

Noting that D is the non-abelian term we see that the non-abelian nature of the field 
is crucial in avoiding the Bekenstein argument for no-hair. The above relation also shows 
that the value of C ■ D must lie in a very restricted range for the non-positivity of S: the 
cubic term in the fields involving the structure constants is bounded both above and below 
by the quadratic and quartic terms in the gauge fields. (It is possible to make restricted 
gauge transformations that preserve time independence of fields and A$ = 0. What we see 
is that (2.5) must hold in some region in each such gauge.) 

The left-hand side of eq.(2.1) also presents complications. For a gauge field, b^b^ is not 
a physical scalar and need not be bounded on the horizon, so a practical application of the 
Bekenstein approach requires the explicit evaluation of 

= --LfW^a^ (2.6) 

on the horizon to make the sign considerations for eq.(2.4) relevant. By introducing the 
ansatz for the connection A used in [19] and [18], it can be shown that b^ is bounded on 
the horizon and that the second term in (2.4) can be negative in the black hole exterior, so 
that a black hole with nontrivial gauge field structure outside the horizon is allowable (but 
not necessary) under the Bekenstein analysis. In fact, the solutions of [18] are examples of 
these "colored" black holes. 

When this approach is applied to the ENAP theory considered below, the results are 
similar. For our metric convention, the Lagrangian is 

2 (i) (i) 

with \A\' = g^A y ;'Ay' and ix is the vector field mass. As we stated above, this is obviously 
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not a gauge-invariant theory but provides a good first approximation for a massive non- 
abelian gauge field. For the massive field, b^b^ is now a physical scalar bounded on the 
horizon, and the Bekenstein integral (2.1) assumes the same form as in (2.4), with £eym 
replaced by £enap- Fot A$ = 0, the same arguments as above establish the possibility of 
hair for massive vector black hole solutions. 



The second theory we investigate in this paper is EYMH, whose Lagrangian is given by 



£eymh = £eym ~~ i — 

4"7T 



[D^ ( W) + V ($)] , 



(2.8) 



where = d^+g r- A M is the usual gauge- covariant derivative expressed in the antihermitian 
representation of su(2) ( n = —iai/2). The Higgs field <3> is taken to be a complex doublet 
$ = (0 + ,0~) and V ($) to be a double-well potential with degenerate vacua. Though we 
give further details of our Higgs ansatz below, for now we assume $ possesses only one degree 
of freedom: 

V2 V<Mx) 

for <fi real and time-independent. The Lagrangian (2.8) then becomes 



(2.9) 



£eYMH = ~t [ \ |F|2 + \ \ 1 



\A\ 2 + l -\d^\ 2 + V^) 



(2.10) 



where \d<$>\ = g^d^dycj) > 0, and (g<j)/2) now plays the role of the gauge field mass. The 
vector b^ has an additional term 



47T 



and the Bekenstein integral is 



(2.H) 



■— / d x^i 



n 



§I F I"+ ( y I W 



9 e ijk-^j) A\j ^ 



lW+(y) W + V'iftt 



(2.12) 



with V' ((f)) = (0) / c?0. The separate integrals in this expression correspond to the YM 
and Higgs surface integrals obtained from eq.(2.11). The YM term is identical in form to the 
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Proca case, so by the above arguments YM hair is possible. In the Higgs term, the sign of 
V (0) for a double-well V (0) obviously depends on the details of a solution, so the second 
integrand need not be positive definite either. Black hole solutions to EYMH theory having 
both nontrivial YM and Higgs field structure are therefore not excluded by the Bekenstein 
analysis. 

It is important to emphasize that these sign arguments only demonstrate the possibility 
of solutions with hair. An example for which all of the above reasoning is valid but a no-hair 
proof has been found [2] is the pure Higgs (Goldstone) theory, which corresponds to the 
Higgs contribution to (2.12) with \A\ = 0. The Higgs-gauge field coupling and the non- 
abelian gauge interactions in EYMH theory therefore play a critical role in the existence of 
black hole solutions for which <p ^ constant outside the horizon. We should also mention 
that Price [4] has employed somewhat different reasoning to argue for a black hole no-hair 
theorem. His arguments, though, are only valid in the domain of weak fields (in which the 
important nonlinear terms would lose their influence) and hence do not directly apply to 
the present context. In the following we shall show that the gap in the no-hair arguments 
described is wide enough to allow for solutions which have spontaneously broken gauge and 
Higgs hair. 

III. ENAP Theory: Regular and Black Hole Solutions 

Metric and Connection 

The metric for a static, spherically symmetric spacetime can be written 

ds 2 = -T- 2 (r) dt 2 + R 2 (r) dr 2 + r 2 {d6 2 + sin 2 9dip 2 ) , (3.1) 

where R (r) = (1 — 2m/r)~^ 2 and m (r) may be interpreted as the total mass-energy within 
the radius r. An alternate form of the metric convenient for describing black hole solutions 
follows from defining 5 = — In (R/T). In terms of the functions (5 (r) , m (r)), the metric 
becomes 

ds 2 = - (l - ^ e- 25 dt 2 + (l " ^) 1 dr2 + r2 K + sin2 * V) . (3.2) 

For the metric (3.1), regularity at the origin requires T (0) < oo and R' (r) ,T' (r) — > , 
while asymptotic flatness requires R(r) ,T (r) — > 1. For the alternate form of the metric, a 
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regular event horizon at r = requires 



m (r h ) = r h /2, 5 (r h ) < oo, (3.3) 

and comparison with the asymptotics of eq.(3.2) gives S (r) — > as r — > oo. For static 
solutions, however, we can rescale the time coordinates 



dt = T- 1 (0) dt 
dt = e~ S ^dt 

in order to simplify the initial data problem: for such a rescaling, 

f (0) = 1, f (oo) = l/T o = l/T(0) 
5(r ft ) = 0, 5(oo) = -(J o = -(J(0). 



(3.4) 



(3.5) 



Thus the correct initial values of gu for a given system can be determined by integrating the 
Einstein equations with the fully specified initial conditions 

i2(0) = l, T(0) = 1 

(3.6) 

m(r h )=r h /2, 5 (r ft ) = 

where in each of the last three equations, the top part refers to regular solutions and the 
bottom part to black hole solutions. We will use this parameterization below. 

The most general spherically symmetric SU(2) connection [26] is 

A— - \a? r dt + br r dr + (d? e - (1 + c) %) dO + ((1 + c) ? e + dr^) sin 6d(p] , (3.7) 
9 

where a and b have dimensions c and d are dimensionless, g is the gauge coupling 

and (T r ,TQ,Tp) is the antihermitian su(2) basis expressed in the usual 3-d (physical space) 
polar coordinate directions; e.g. % = r • r, and [r a , %] = t a bt^c with the indices ranging over 
(r, 9, <f). The four degrees of freedom a, b, c, d are all functions of the 3-d radius r and time 
t. The connection (3.7) has a residual gauge freedom 



UAl/- 1 + -Udl/- 1 (3.8) 



under unitary transformations of the form U = exp [(3 (r, t) %], where (3 (r, t) is an arbitrary 
real function. Under such gauge transformations, which form an abelian subgroup of the full 
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gauge group, the connection functions transform as 



U : 



b 

c 

\dj 



b 

c 



( a-P \ 
b-(3' 
c cos (3 — d sin (3 
\ d cos (3 + c sin (3 ) 



(3.9) 



We can use this freedom to impose the "polar gauge" 6 = 0. In the static case, (3 = and 
a, c, d are functions of r only. Following [19] and [18], we eliminate two of the three remaining 
degrees of freedom by the ansatza = (no dyons, purely magnetic YM curvature) and d = 0. 

After relabelling the remaining degree of freedom w (r) = c (r), the connection assumes 
the form 



A = - (1 + w) [—TtpdO + rg sin 9dip] 



(3.10) 



which was explored by t'Hooft in the context of magnetic monopoles [28]. It differs from the 
ansatz of [19] by a singular gauge transformation U = exp {9t\) exp ((f — V 9 ) r 3) an d has the 

2 

added virtue that |A| is invariant under spatial rotations. 
ENAP Equations and Boundary Conditions 

Defining F = dA + gA A A , we have 



w 



1 



F = — \-%dr Ad9 + ?g sin A dtp] (l - w ) % sin 6d6 A d<z>. 

9 9 

The Proca equations D*F + fi 2 *A = reduce to a single equation: 



(3.11) 



/ w (1 - w 2 ) i? 2 

rfr V RT J 



2 y-^a+^f = o. 



(3.12) 



The Einstein equations = ftnT^p may be calculated by varying with respect to g^j, 
the action 

J d 4 x^ (JL + £ ENAP ) = ^- J d 4 x^ (R - \F\ 2 - 2/i 2 |A| 2 ) , (3.13) 
where R = R iiv g {iv is the gravitational curvature scalar and we have set G = 1. The 
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energy-momentum tensor is 



87rT M „ = 2F^FS - X -g» v |F| 2 + 



2 l ?A tl A v -g llvl ?\A\ 2 \ , (3.14) 



and the (tt) and (rr) Einstein equations can be written 



m 



2m\ fw'V l (l-w 2 Y u 2 
r )\g) + 2 g 2 r 2 + g* ^ 

f 2m\ fv/\ 2 1 jl-w 2 ) 2 
V r J \g J + 2 ^ r 2 



+ w) 



(3.15) 
(3.16) 



+ ^(i + w y — 



in terms of the metric functions m (r) and T (r). While (3.16) is appropriate for obtaining 
regular solutions, the coordinate singularity at the event horizon of a black hole solution 
leads us to replace it with 

.2 

(3.17) 



8' 



2 (w'Y 



9 ' 



where 5 was introduced in (3.2) above. For either choice of metric functions, the two Einstein 
equations can be used to express the Proca equation (3.12) in a form independent of T (r) 
or 5 (r): 



2 / , -"'\ // 
r I 1 w 

r J 



2m- 



-2%(l + w) 2 r 



g 2 r 



9 



w 



(3.18) 



+ (1 - w 2 ) w - jj 2 (1 + w) r 2 = 0. 



All of the above agrees with the corresponding EYM results [19,18] when we take \i = 0. 

The presence of the additional parameter // in ENAP theory motivates us to clarify the 
role of the gauge coupling g in the field equations. A classical theory with G — c — 1 
satisfies [L] = [T] = [M], which for our theory implies [g] = [/j] = IX] 1 . If we introduce the 
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dimensionless quantities 



r = gr, m (?) = g m (gr) , Jl = fi/g, (3.19) 

then the theory scales as 

\F\ 2 = g 2 \F(r)\ 2 , \A\ 2 = \A(r)\ 2 , £ E NAP = <? 2 £enap (r, $) , (3.20) 

and the dimensionless ENAP equations assume the form of (3.15)-(3.18) with g = 1 and 
(3.19) replacing the dimensionful variables and parameters. Through scaling, we may there- 
fore obtain a solution to the ENAP equations for any g > from a solution to the dimen- 
sionless equations: 

w g (r) = W (gr) R g (r) = (1 - 2m (gr) / (gr)y 1/2 

1 (3.21) 
m g (r) = -m (gr) T g (r) = T (gr) ; S g (r) = 5 (gr) . 

Solutions with g > 1, for example, possess the same structure as g — 1 solutions, but it 
occurs at radius r = r/ g < r. For the remainder of this section, we take g — 1 without loss 
of generality. 

We can anticipate from the field equations the general characteristics of solutions. From 
the definition of 5, eq.(3.17) demonstrates that R/T increases monotonically with radius. T 
can also be shown to satisfy the same equation as in EYM theory [19], 



d_ 

dr 



2 / i \ > 



R \T 



2(l-^) M 2 f+(^)(l-^, (3.22) 



so that T' < for w ^ 1. From the boundary conditions of the unrescaled metric (T (0) > 
R(0) = 1, T (oo) = R(oo) = 1) and the definition of R, we then expect that T > R > 1, 
R should possess at least one maximum, and T should decrease monotonically for nontrivial 
regular solutions. For black hole solutions, 5 should decrease monotonically and R should 
decrease from its singular value at the horizon, though it need not possess any local extrema 
for r/j < r < oo. 
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Following [19], we can learn about the behavior of the vector field by rewriting the Proca 
equation in the form 



1 d 



RT 



+ ^- [(w 2 - l) w 2 + fi 2 (I + w) wr 2 ] . (3.23) 



RT r 2 T 



A trivial m = solution occurs for w = — 1, and comparison with eq.(3.15) reveals that 
w = — 1 is the only acceptable asymptotic value for finite energy solutions. Since the right- 
hand side of this equation is manifestly positive for w 2 > 1, the only nontrivial solutions 
having finite w 2 and finite energy must satisfy w 2 < 1. The precise finite energy restrictions 
on the behavior of w come from the nonderivative term in (3.23), and since w' — gives the 
relation 

R 2 

r 2 



ivw" = ^ (1 + w) w [w (w - 1) + /i 2 r 2 ] (w' = 0) , (3.24) 



this term also governs the oscillatory properties of finite energy solutions. From (3.24), 
we see that w" > for r > 1/2//, so finite energy solutions must have their final turning 
points before r = l/2/i and satisfy w' < after r = 1/2/x. For r < l/2/i, on the other 
hand, eqn.(3.24) indicates that ww" < in the range ^1 — yj 1 — (2yur) 2 ^ /2 < w < 

1 + y 1 — (2/ir) 2 ) /2, so solutions can exhibit nontrivial behavior for w > and still 



satisfy the boundary conditions at infinity. If w > ^1 + 1 — (2fir) 2 ^j /2, however, w will 

increase beyond u> 2 = 1, so this range is forbidden for finite energy solutions. Since we have 
used only the r — > oo boundary conditions and not mentioned initial conditions on w or the 
metric in this discussion, the results 



r<i-: -l< w; <i^l + v /l-(2 / ir) 

1 / // 

r > — : w < 0, «/ > 
2/i 



(3.25) 



describe both regular and black hole solutions. 

The full boundary conditions may be determined from eqs.(3.15)-(3.18). For regular 
solutions, finite energy density Tu and regularity of T (r) at the origin (T' (0) = 0) give 
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2m (r) = O (r 3 ) (3.26) 
lnT(r) = 0(r 2 ) (3.27) 
w (r) = -1 + (r 2 ) . (3.28) 

It is significant that w — — 1 is the only possible initial value; in EYM theory the self- 
interaction term (l — w 2 ) jr 2 in w! provides the two possibilities w± = ±1 [19], but the 
vector field mass in our theory excludes w + both at the origin and as r — > oo. The fj 2 
terms also breaks the discrete symmetry w — > — w of the EYM equations which gives rise to 
degenerate mirror-image solutions in [19] and [18]. These differences will play an important 
role in the spectrum of solutions. 

For black hole solutions, the field equations with m (77J — r^/2 give m! (rh) , w' (r/J, and 
8' ij'h) when jj, is specified. Full use of the metric initial conditions allows us to expand near 
the horizon: 

m (r) = r h /2 + m! (r h ) (r - r h ) + O (r - r h f (3.29) 
5 (r) = + 5' (r h ) (r - r h ) + O (r - r h f (3.30) 
w (r) = w (r h ) + w (r h ) (r - r h ) + O (r - r h ) , (3.31) 

where and w (r^) are to be chosen to yield a finite energy solution. 

Since the only allowable vacuum value of the connection function for both regular and 
black hole solutions is w = —1, the behavior of eqs.(3.15)-(3.18) as r — > 00 gives 

m (r) ~ M - fic 2 e- 2fir (3.32) 
InT(r)~ln^+^ (3.33) 

5 (r) ~ -6 + fic 2 , (3.34) 

r 

w (r) 1 + ce^ r , (3.35) 

where c is some positive constant and T and 8 refer to the unrescaled initial values intro- 
duced in eqs.(3.5). Thus the presence of the fu 2 terms in ENAP theory give exponentially, 
rather than polynomially, decaying fields as we expect. 
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Numerical Regular Solutions 

We use a standard one-parameter "shooting" method to find regular solutions to the 
ENAP equations. The formal power series describing the boundary conditions (3.26)- 
(3.28) at r = is 

2m (r) = 4&V + - (-8b 3 + 3b 2 /i 2 ) r 5 + O (r 7 ) (3.36) 
InT (r) = -2&V - - (12& 4 - 46 3 + fi 2 b 2 ) r A + O (r 6 ) (3.37) 
w (r) = -1 + br 2 + ^ (86 3 - 3b 2 + /x 2 6) r 4 + O (r 6 ) , (3.38) 

which depend only on b > and and conform to the EYM results [19] for fi = 0. We 
use these conditions evaluated at r = 1CP 3 as initial data in a standard ordinary differential 
equation solver with a global error tolerance of 10 -12 , adjusting b for fixed jj, and integrating 
outward in an attempt to meet the boundary conditions at infinity. The bracketing condition 
for finite energy solutions is the same as in [19]: for a small range of b in the vicinity of a 
solution, w approaches its asymptotic value w — > — 1 and either crosses through w — — 1 
and rapidly goes to w = — oo, or experiences a turning point and rapidly goes to w — +oo. 
Though 2m (r) — > r as diverges, at a discrete value of b in this range lies a finite 
energy regular solution with the correct asymptotics. The existence of such solutions has 
been proven rigorously only in the \x = case [22], but the qualitatively similar numerical 
behavior for ji > lends strong support for our work. 

As in EYM theory, finite energy solutions are characterized by the oscillatory behavior 
anticipated above in the near-field region r ^ 1. Distinct solutions are classified by k, the 
number of zeros of w, which increases with increasing b but must be even for \i > to 
conform to the boundary conditions (3.35) and (3.38). A striking difference between our 
non-abelian Proca solutions and the EYM solutions of [19] is that for /i > 0, we find two 
distinct classes of even-A; solutions. We will now discuss each class in turn. 

One class possesses physical characteristics very similar to the even-A; solutions of [19] 
for small //, reducing to them in the limit ji — > 0; the presence of the vector field mass only 
slightly perturbs the gross behavior of the field. Because it offers the best hope for dynamical 
stability, we focus our attention on the k — 2 solutions. In the range < fi < 4.454x10^2, the 
shooting parameter varies over 0.6517 > b > 0.5787 while the width of the single peak of w (r) 
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increases as the mass increases: 0.9713 < M < 0.9936 (cf. fig. la). For \i > 4.454 x 10 , 
w does not approach its asymptotic value before an additional turning point occurs and 
w — > +oo: the /i 2 term in the Proca equation begins to dominate at a radius r ~ l/fi fa 20 
before the non-abelian self-interaction becomes negligible, thereby making impossible A;-node 
solutions for any value of b. Another way to understand the occurrence of a maximum ji 
value is through eqn.(3.15); if the /i 2 term dominates the energy density, the radius will 
always approach the Schwarzschild radius r = 2m beyond some value of /i and solutions 
with smooth geometry will be forbidden. 

The second class of solutions resembles the (k — 1) odd-node solutions of [19] as /t — > 0, 
with w approaching the forbidden asymptotic value w — +1 until large r < 1/2//, where 
a turning point and then and additional node occur before w — > — 1 (fig. lb). As the 
mass parameter increases toward \x = 4.454 x 10~ 2 , this type of solution begins to resemble 
the limiting case of the first solution type, until eventually the two classes converge at the 
maximum value of /t: for < /t < 4.454 x 10~ 2 , the shooting parameter and total mass 
vary over 0.4537 < b < 0.5787 and 0.8286 < M < 0.9936, respectively. Thus solutions to 
ENAP theory bifurcate into two branches, with the bifurcation point corresponding to the 
H at which /c-noded solutions no longer occur. In shooting parameter space, the bifurcation 
condition corresponds to the shrinking of the interval over which k nodes occurs from the 
maximal value of [19] (0.4537 < b < 0.6517 for k — 2) to zero, with finite-energy solutions 
occurring at the endpoints. 

The existence of a second class of solutions and the bifurcation phenomenon can be 
understood from a heuristic argument of scales. In EYM theory, there is one dimensionful 
parameter g whose units are [T] [M]" 1 / 2 [L] -3 / 2 (where for the remainder of this paragraph 
we do not set G = 1) which sets the scale of the solutions to be on the order of Gr 1 / 2 ^ -1 . 
By scale here we refer to the approximate value of the radius beyond which the fields in 
the theory exhibit no nontrivial behavior and rapidly approach their asymptotic values. In 
ENAP theory, on the other hand, we have two dimensionful parameters [g] and \p] , the latter 
of whose units is [L]^ 1 . There are thus two distinct length scales : G 1 / 2 ^ -1 and fi^ 1 . The 
former sets the scale for one class of our solutions and the latter sets the scale for our second 
class. Notice that the first scale is dependent on gravitational interactions while the second 
is not. As I /n — > oo, the last node of a A;-node solution in the second family is pushed 
off to infinity and the solution approaches an odd-node solution of the EYM system with 
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its scale being set by the only remaining parameter, G 1 / 2 g~ 1 . These quasi-odd-A; solutions 
[29] therefore have the scale of their inner structure set by gravity while the scale of their 
asymptotic structure is set by the Proca mass [31]. 

Numerical Black Hole Solutions 

To find numerical black hole solutions, we followed a similar shooting procedure. The 
field equations give 

rri (r h ) = [l-w 2 (r h )] 2 /2r 2 h + ^ [1 + w (r h )f 
w , ^ _ y 2 [l + w (r h )} r 2 h -[l- w 2 (r h )] w (r h ) 

r h - [1 - w 2 (r h )} 2 /r h - 2^ 2 [1 + w (r h )} 2 r h 
5'(r h ) = -2[w'(r h )] 2 /r h 

on the horizon, and we can construct initial data from eqs.(3.29)-(3.30) with only the shoot- 
ing parameter w (r^) and the horizon radius r\ unspecified. Following [18], we examine 
rh = 1 and chose (r — r^) < 1CT 12 so that the errors in the initial data are smaller than the 
global tolerance. 

As in the smooth case, we again find two distinct classes of solutions whose scales are 
respectively set by the two dimensionful parameters in the theory. In appropriate limits, these 
solutions approach those of [18] modulo one interesting special case to be mentioned below. 
For n > we find that odd-node solutions are possible in addition to even-/c solutions, since 
the horizon shields the singularity which would occur at w (0) = +1 for such regular solutions. 
These /c-node solutions are classifiable by their behavior as /i — > 0; one class reduces to the 
A;- node solutions of [18], while the other approaches the (k — l)-node solutions of [18] as the 
position of the /cth node moves out to r = oo. The two solution branches are again joined 
at a bifurcation point for some maximum value of fi for each k. The k = 1 quasi-even-node 
case is special because there exists no corresponding non-abelian solution with A; — 1=0 
nodes in [18]; instead, this solution approaches w = +1 as fj, — > 0, which corresponds to the 
ordinary Schwarzschild solution. The (regular and black hole) analogue of this solution in 
Einstein- Yang-Mills-Higgs theory (to be discussed shortly) is significant as it allows us to 
make contact with the known existence of a flat space sphaleron. This is the only set of 
solutions (black hole or regular) for which the limiting case is essentially weak-gravity (weak 
in the present sense meaning that it becomes arbitrarily close to Schwarzschild). Hence, 
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(3.39) 
(3.40) 

(3.41) 



ENAP black holes are a set of fundamentally non-abelian solutions which, like the solutions 
of [18], possess metrics which interpolate between the RN and Schwarzschild metrics, but 
which include as limiting cases both quasi-Schwarzschild and quasi-RN (k — > oo here and 
in [18]) solutions. In this sense, the spectrum of solutions itself interpolates between the 
Schwarzschild and pure-magnetic RN black holes. 

Though the solutions for all node numbers exhibit bifurcation, we focus our attention 
on the lowest odd- and even-node solutions. In figs. 2a-b, the mass and connections for 
the two branches of the k = 1 solutions are plotted. The limiting value of fi = 0.1233 
joins the solution branches, which are described by 1 > w (rh) > 0.8500, 0.5000 < M < 
1.0052, < 5 < 0.2854 (quasi-A; = 0) and 0.6322 < w{r h ) < 0.8500, 0.9372 < M < 
1.0052, 0.5485 > 5 > 0.2854. For k = 2, the maximum value fx = 1.752 x 10" 2 gives the 
branches -0.6322 < w (r h ) < -0.5027, 0.9372 < M < 1.0052, 0.5485 < 5 Q < 0.5831 (quasi- 
k = 1) and -0.3425 > w (r h ) > -0.5027, 0.9938 < M < 1.0052, 0.5932 > 5 Q > 0.5831 (cf. 
figs. 2c-d). For small /i, the deviation from the Schwarzschild mass M = 0.5 appears to 
scale with /i 2 for the quasi-A; = solutions; this can be understood from (3.39) , which for 
w ^ +1 gives m! (r/>) ps 4/i 2 independent of r^. 

As we vary — > 0, we find that the even- A; black hole solutions reduce to the regular 
solutions for the same value of /i in the manner of [18], with the shooting parameter ap- 
proaching w (r/J = — 1 +b r^, where b is the regular solution shooting parameter. The value 
of ji at which the two solution classes bifurcate also increases as — > from its value at 
rh — 1 to the regular solution bifurcation value; for k = 2, 1.752 x 10~ 2 < fibif < 4.454 x 10~ 2 
as < r/j < 1. The reduction to the regular solutions is sensible when we examine eqs.(3.39)- 
(3.41); despite the presence of the horizon, w' (rh) and 8' (rh) approach the leading behavior 
of the r — > regular solution expansions with r replaced by rh, while the non-horizon mass 
contribution (2m (r) — rh) /(r — rh) = 4 6 2 r^ mimics the regular solution mass-radius ratio 
2m (r) jr in the same limit. 

Although for fj, ^ there are no regular solution limits which the odd-A; black holes can 
approach as r/j — > 0, there is limiting behavior as the horizon continues to shield the w (0) = 
+ 1 singularity. The shooting parameter for such solutions behaves as w (rh) = +1 — b e q r^, 
where 6 e ff is now an effective regular shooting parameter, and we can determine bifurcation fi 
values as well as mass limits for the solution branches. In the k = 1 case with r/j = 10~ 3 , for 
example, we find for < \i < 0.2010 the branches < b eS < 0.1329, < M < 0.9289 (quasi- 



19 



k = 0) and 0.4537 > b cS > 0.1329, 0.8286 < M < 0.9289; the latter branch reduces to the 
k = 1 solution of [19] while former approaches the Schwarzschild solution with M decreasing 
proportional to r\. Because the curvature at the horizon = SirT^ = /i 2 (l +w) 2 /rf i 
diverges as — * 0, the maximum value of /x will presumably decrease until it falls off 
faster than ~ r^, which will bring the expansions (3.40)-(3.41) and the non-horizon mass 
contribution into the regular solution form (we have not investigated this in detail). Thus 
the only legitimate quasi-regular limiting black hole solutions for odd k occur for fj, — > 0, but 
these are significant: the — > odd-node black holes help "complete" the regular spectrum, 
in that the regular and quasi-regular solutions to ENAP theory interpolate between the 
Schwarzschild vacuum and the gravitating Dirac monopole [17]. As mentioned above, we will 
find that the quasi-/c = branch of solutions, which are trivial apart from a rapid transition 
between the minima of the non-abelian self-coupling w — +1 and w — — 1, exist with or 
without event horizons in EYMH theory and correspond to gravitating SU (2) sphalerons. 

It should be noted that the divergent behavior of w as r — > oo used to shoot both 
regular and black hole solutions is characteristic of integration toward a singular point. We 
would observe the same qualitative behavior as r — > or r\ if we had chosen to integrate 
inward from some large value of r using (3.32)-(3.34) as initial data and {c, M} as shooting 
parameters. To help verify the existence of our solutions, we integrated from both small r 
and large r to a point inbetween, where we attempted to match the values of w, w', and 
m by adjusting the shooting parameters {b,c,M}. Using the value of b found from our 
original method, we were able to adjust {c, M} to give agreement at the common point to 
an accuracy comparable to the global tolerance. We repeated this procedure for EYMH 
solutions after including additional parameters associated with the Higgs field. 

IV. EYMH Theory: Regular and Black Hole Solutions 

The Higgs Ansatz 

The starting point for our theory is the addition in eq.(2.8) of a complex scalar La- 
grangian to the ordinary YM Lagrangian. Following the standard model, we take the Higgs 
field to be a complex doublet $ = (0 + ,</>~) with the double- well self- interaction 

V($) = \ (V$) - /i 2 + constant. (4.1) 
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The most general complex doublet Higgs may be parameterized 



1 / 

$ 0*0 = -7= exp [-r • £ (x)} 

V2 \h(x)/r 



(4.2) 



with /i and £ a arbitrary. A parameterization of the Higgs field which proves useful in finding 
spherically symmetric EYMH solutions can be achieved by rewriting $ in the form 



x 



l (ih(x) + #i 

$(x) = —=\ , . 

V2 y - #3 (a) / 



(4.3) 



where we have grouped three of the degrees of freedom as a vector ip. Substituting into the 
EYMH Lagrangian (2.8) gives 



£eymh = - 



1 

47T 



l lTll2 1 

4 1 1 8 



+ 



+ 2<T + (^) • + ^ (V + l^l 2 ) ( 4 - 4 ) 



where A is our original connection ansatz (3.7) . A useful ansatz for the Higgs field which 
yields a spherically symmetric energy density is to take <fi = <fi (r, t) and if) — ip (r, t) n r , with 
n r a unit vector in the radial direction. 

The static field equations in the gauge 6 = 0, without the ansatz a = d = 0, are 
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[1 (^ + ^)] ^ _ 1 [(1 + w) 2 _ (1 _ w) ^ p + fl 2 wi?T = (46) 

(4.7) 



T 4 

[l — (t(7 2 + G? 2 )] i? 1 



d- - - [(0 2 + ^ 2 ) d - 2^} - + a 2 dRT = 



vy- - - 

<t> T 2 
,i? 1 



1 + w f + d 2 - 2#) - + -a 2 r 2 (t)RT = (4.8) 
^> 2 ^ - \ ( [(1 - ™f + d 2 ] V> - 2rf0) | + iaVtyiZT = 0, (4.9) 
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with the constraint equation 



2 
RT 



(d'w - w'd) + -r z 
2 







(4.10) 



arising from the gauge choice. Following Witten [26], it is useful to express the w and d 
degrees of freedom in complex scalar form: 



w (r) — id (r) = / (r) exp ia (r) . 



(4.11) 



Substitution of w = f cos a and d = —f sin a into the constraint equation with = ^ = 0, 
v = gives the result / 2 c/ = 0, implying that the gauge choice allows for the dynamical 
elimination of one of the remaining degrees of freedom: w and d are related by a multiplicative 
constant (the equations for d and w in EYM theory are identical), and d = in [19] and 
[18] is not an ansatz but a result of the choice of the constant a = [19, 18]. With and ^ 
nonzero, the constraint equation becomes 



2/V + i ^ 7 



2./ 



0. 



(4.12) 



Examination of the field equations then implies that if either a or 7 is chosen to be a constant, 
then a — 7 = nn. 

We will only study the case in which a and 7 are constants, and hence obtaining static 
solutions to EYMH theory reduces to solving the coupled equations for a,f, and h/r and 
choosing a and n. 

The equation for a may be rewritten 



1 d 

2dr 



2^ I 2\' 



,T 



— (a) +2RT 
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while the / and h/r equations become 
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(4.13) 



r [f + co^Y- + a 2 fRT 



R 1 



h\ R 1 
~r T + ~^ r 



2^2 



h 



RT = 0. 



(4.14) 
(4.15) 



For odd n, the above equations can be put in the same form as the even-n equations by 
defining / = — /, so the solutions for odd and even n are isomorphic despite the fact that 
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the physical field configurations are quite different. Before exploring these differences, we 
note that the r — > oo boundary conditions on / and h require a (oo) = 0, and for regular 
solutions a (0) = 0. Since the right side of (4.13) is positive semi-definite, we must have 
a = for static regular solutions. Thus the no-dyon results for EYM theory [32] extend to 
EYMH theory. 

Though solutions to (4.14)-(4.15) do not depend explicitly on the values of a and n, 
some physical properties of the Higgs field and connection are affected by these parameters. 
Consider a — 0, for example. The even-n solutions reduce to our ansatz, while the odd-n 
solutions have have a "hedgehog" Higgs field ((f), if)) = (0, ±/i/r). We will shortly see that 
one of the smooth k — 1 node solutions of this form which we find in the next section has a 
limit in which it becomes close to the flat space sphaleron of YMH theory. 

EYMH Equations and Boundary Conditions 

To obtain static, spherically symmetric solutions we assume = (r) = v + rj (r) and 
■if) = 0. The Lagrangian for the theory may then be written in the form 

£eymh = ~ (l W\ 2 + \ ( y) Vl 2 + \ \H 2 + V (0) j , (4.16) 

where V ((f)) = A [(f) 2 - v 2 ] 2 /4 satisfies V (±v) = 0. 

Specializing our results of the last section, the gauge field equation is 

dfw'\ w(l-w 2 )R (g(f)\ 2 ,, ,R n , Ain , 

and the Higgs equation is 

±(*t-)-l(l + w) 2 <f>*- V ((f)) r 2 * = 0, (4.18) 
dr \RT J 2 ! V T T K ' 

where V = dV/d(f) = dV/drj. 

From variation of the full action with respect to we obtain the energy-momentum 
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tensor, 



8tt7V = 2F^FS~g^ \F\ 2 + 



+2d^d l/ r ] -g lll/ \dr 1 \ z + 2V{<P) 



(4.19) 



and the (tt) and (rr) Einstein equations: 
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(4.21) 



For black hole solutions, we again replace the auxiliary T' equation with an equation for 8', 



5' 



2 2~ 

2 (w'/p) + (W) l r - 



(4.22) 



The Einstein equations can once again be used to express the equations of motion in a form 
independent of gtt' 



(l-w 2 ) 

2m -± T-^~ 2 [-\ [l + wY r - 2V (<p) r £ 

g r \ 2 
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(l-w 2 )w-(^p) (l+w)r 2 = 



(4.23) 



r 2 [l- 2 -^)v" + 



2 (r — m) — 



(1-w 2 ) 



g*r 



2i -2 ( I ] (l + u>) r-2\/(0)r ; 



- (1 + wf<j)- V'{<t>) r 2 = 0. 



(4.24) 



The addition of the Higgs field and new coupling constants A, jj 2 requires a reexamination 
of the scaling properties of the field equations. From the replacement of /i 2 in the ENAP 
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theory with (g<f>/2) and examination of £eymh> we find 



fo] = M = l, [A] 1 / 2 = M = [L]- 1 , 



1-1 



(4.25) 



so by introducing X = X/g 2 and /2 = /x/gr, we obtain the same overall scaling behavior of the 
Lagrangian found in section III: 



£eYMH = fi^EYMH (r, A, fij 



(4.26) 



Solutions for g ^ 1 will again be related to dimensionless quantities as before, with the 
addition of r\ g = n (gr), and the dimensionless field equations can be obtained by replacing 
(r,m,X,/i) with (r, in, A,/2^ and setting g — 1 in eqs.(4.20)-(4.24). We take g = 1 for the 
remainder of the paper and consider all quantities dimensionless unless otherwise specified. 

By following the analysis of section III, we can again predict general features of the 
solutions and boundary conditions. From the expression for S (4.22) , we see that R/T again 
increases (and S for black holes decreases) monotonically with radius, but because T now 
satisfies 



d_ 

dr 



^©l-( 1 -v)M 2 f + (^)[( 1 -^-- W ], (4,7) 



it is no longer clear whether T decreases monotonically for regular solutions. The condition 
T > R > 1, with R possessing at least one maximum, still holds for regular solutions. 

The YM equation may be rewritten 
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(4.28) 



in general, while at turning points we have the relation 



R 



ww" = (1 + w) w 



w (w - 1) + ( - I i- 



(w' = 0) 



(4.29) 



Since <fi (oo) ^ for a spontaneously broken theory, w = — 1 is again the only acceptable 
asymptotic value and w 2 < 1 is required for finite energy solutions. Because is a field 
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and not a constant, determining restrictions on the value of w and the occurrence of turning 
points is not as simple as for ENAP theory. If we assume that is O (v) in the region 
r <; 1, we can again use a characteristic radius to define two regimes of interest and apply 
the analysis of section III: 



r<0{l/v): -Kw<i^l + y / l - (vrf 
r >0(l/v) : w' < 0, w" > 



(4.30) 



quantify the finite-energy constraints on both regular and black hole solutions. 

To better understand the expected behavior of 0, we rewrite the Higgs equation 
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(4.31) 



The obvious requirement for finite energy solutions is V' (0) < 0, which implies that is 
restricted to lie between the minima of the potential: |0| < v. We also note that = ±v , 
are the only allowed values as r — > oo, but because V (0) is nonzero at (oo) = 0, finite 
energy solutions must have (oo) = ±v. From (4.31), the equation governing the oscillatory 
properties of is 



00" = R 2 



(1 + wY 
2r 2 



2 + V (0) 



(0' = o) 



(4.32) 



When the gauge field coupling term is negligible, we see that 00" < 0, which is characteristic 
of oscillations about = (a solution of infinite energy) unless the initial value of <$' in 
(4.31) is large enough for the field to reach one of the minima of V (0) as 0' — > 0. For 
regular solutions, we will find that 0' (0) = for (0) 7^ 0, so oscillatory behavior will occur 
unless the gauge field coupling becomes important. For black hole solutions and regular 
solutions with (0) = 0, the initial derivative of is nonzero, but the gauge field term again 
plays an important role in avoiding infinite energy solutions. The behavior of the scalar 
field is consistent with the analogue of a particle moving in a potential proportional to —V, 
with the gauge field interaction and gravitational forces fighting the restoring force of the 
potential. These heuristic considerations are in agreement with the no-hair proof for the 
Goldstone theory [2], which approaches the existence of finite energy solutions rigorously 
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but for which the analysis fails in the presence of a non-abelian gauge field. Thus we expect 
the Higgs field either to smoothly transit the potential until \<f>\ = v, or to move initially 
toward = and then turn and roll under the influence of the gauge field toward \(f>\ = v. 
The — * —0 symmetry of the field equations also implies that every solution will have a 
mirror-image solution; we focus here on solutions with (oo) = +v without loss of generality. 

The presence of the Higgs field in EYMH theory allows for two possible sets of initial 
conditions for regular solutions. Finite Tu and regularity of the metric at r = give 



2m (r) = O (r 3 ) (4.33) 
lnT(r) = 0(r 2 ) (4.34) 

Even - k (4.35) 
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Odd - k, (4.36) 



where k again denotes the number of w nodes and — v < i] < for (oo) = +v. Black 
hole solutions again possess expansions like eqs.(3.29) - (3.31) near the horizon, with the 
addition of 

V{r)=v (r h ) + V' (r h ) (r - r h ) + O (r - r h f , (4.37) 

where r\ (rh) brings to three the number of unknowns. 

The vacuum values w (oo) = — 1 and rj (oo) = are shared by black hole and regular 
solutions; the behavior of the field equations as r — > oo gives 



m (r) 
In T(r) 

S(r] 
w (r] 



M 



l -a 2 (V2^ 



re 



-2V2fir 



r 



To, 

„_l + ce -(«/2)r 



X , 1 2 

~ -d + -a 



~ ae 



(4.38) 
(4.39) 

(4.40) 
(4.41) 
(4.42) 



27 



where c and a are constants, v/2 is the gauge field mass after symmetry breaking and y/2fi 
is the Higgs field mass. 

Numerical Regular Solutions 

The presence of the Higgs field makes the solution of the EYMH equations a two- 
parameter shooting problem. For even-/c solutions, the boundary conditions (4.33)-(4.35) be- 



come 



2m (r) = (Ab z + -V )r 6 + 
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where V , V' Q and V^' are the potential and its derivatives with respect to <fi at <fi = <p = f +^ , 
and the shooting parameters are 6 > and r/ c . For odd-A; solutions, the expansions are 



2m(r)= Ub 2 + ^V + ( V ' ) 2 ^ r 3 + 
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r 5 + C (r 7 ) 

(4.48) 
r A + (r 6 ) 



u, ( r ) = +1 - br 2 - i_ - 36 2 + 4F 6 - X - (r]' ) 2 [1 - 86]^ r 4 + C (r 6 ) (4.49) 
»7 (r) = i/ e r + (s& 2 -26 + 3 (^) 2 + +?V + ^ r 3 + (r 5 ) , (4.50) 

where 77^ replaces r\ as the second shooting parameter. We again evaluate the initial con- 
ditions at r = 10~ 3 and use global error tolerance 10~ 12 in a standard ordinary differential 



28 



equation solver, adjusting either (b,i] ) or (b,rj' ) for fixed A and /i and integrating toward 
r = oo. The finite energy solution bracketing condition for b is similar to that of ENAP 
theory: for r] or rf in the vicinity of a solution, there exists a range of b for which w — > — 1 
but then rapidly approaches |w| = oo, with a discrete value of b giving the correct asymp- 
totics for w. As we adjust rj or r]' , the Higgs field either passes through rj = and diverges 
or experiences a turning point and begins oscillating about r\ = —v, with discrete values 
of (o, i] ) or {b, r)' ) giving the monotonic approach of both fields toward the boundary val- 
ues as r — > oo. The difficulty of the two-dimensional shooting problem was compounded 
by the presence of the two free parameters A and /i. Examination of (4.41)-(4.42)reveals 
that a significant disparity between \/2/j> and v/2 makes the accurate determination of solu- 
tion shooting parameters very difficult, since both fields must simultaneously approach their 
asymptotic values for both bracketing conditions to occur. As a consequence, we focus on 
solutions with \[2[i = v/2 (corresponding to A = 1/8). 

The general properties of finite energy regular solutions are the same as for ENAP theory. 
Solutions are again characterized by w oscillations in the region r ^ 1 and classified by the 
node number k, which may now be odd as well as even. For odd k, the Higgs deviation 
monotonically increases in the range — v < i] < 0, while for even k it undergoes a turning 
point very close to r\ and then monotonically increases to r\ (oo) = 0. Two distinct sheets 
of solutions for each k again arise (for any choice of A or jj, in the appropriate range we have 
two solutions for each choice of k), but for A = 1/8 they do not precisely converge, and 
a different maximum value fmax 

for each sheet may be determined. Though we have not 
fully investigated this phenomenon, these two sheets corresponding to the different solution 
classes presumably bifurcate at a point with A^j ^ 1/8. The results for k = 1 and k = 2 are 
shown in fig. 3a — d. As v varies in the range < v < 0.599 for the first k = 1 branch, the 
shooting parameters vary over 0.454 > b > 0.261 and 0.323t> < rj < 0.329t> max and the mass 
behaves as 0.8286 < M < 0.9272; the v -> limit is the k = 1 solution of [19]. The quasi- 
k = branch is described by < v < 0.619, 0.292w 2 < 6 < 0.206, 0.434w 2 < rj' < 0.307w max , 
and 1.821f < M < 0.9380. This approaches the weak-gravity regime quickly as v decreases: 
gu (0) increases from its minimum value 0.1557 at v umx to g tt (0) = 0.9998 at v = 10~ 2 . We 
also note that the scaling of the mass with v is characteristic of sphaleron solutions to YMH 
theory, but address the relationship between our solutions and sphalerons more fully below. 
For k = 2, the two solution branches are described by < v < 0.120, 0.652 > 6 > 0.585, 
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-.870v < T) < -.845v max , 0.9713 < M < 0.9972 and < v < 0.122, 0.454 < b < 0.568, 
-v < r} < -.850w max , 0.8286 < M < 0.9982 (quasi-A; = 1). As v -> 0, the branches 
approach the k = 2 and k = 1 solutions of [19], with the initial Higgs deviation of the 
quasi-A; = 1 branch approaching the odd- A; initial value r\ — — v. 

We can again understand in heuristic terms the bifurcation of the solution space into 
two branches by appealing to length scales [33] . The coupling parameters in this model 
provide us with two length scales: L\ = g~ 1 G 1 / 2 and L2 = g~ 1 ^/\/ fi 2 (where we do not 
set G = 1 for the remainder of this paragraph). As in our discussion of ENAP theory, the 
former of these scales is set by gravity while the latter is not. One branch of our solutions 
has its scale (the position after which the fields rapidly approach their asymptotic values) 
set by L\ and the latter by L2. In the limit that the v goes to zero the k node solutions on 
the latter branch have their last node pushed out to infinity and they go over to k — 1 node 
solutions with scale set by L\. An interesting special case is that of our k = 1 node solutions 
(smooth or black hole) since they would appear to approach a k = node solution. This 
means that their scale is set totally by L2 and hence can exist with arbitrarily weak gravity. 
In particular, our smooth solution of this sort can be related to the flat space sphaleron in 
this limit. Explicitly, taking our solution with the choice a = and odd-n we can write the 
solution in the form 



V2\ 



(4.51) 



U = exp [—miT r } , 



which is physically equivalent to the YMH sphaleron ansatz of [34] and [35]. Though the 
sphaleron solutions utilizing (4.51) have only been found for k — 1 in YMH theory, we see 
that they appear to exist for all k when we include gravity. 

Numerical Black Hole Solutions 

To find numerical black hole solutions, we used the conditions 

2 

[l+ W (r,)] 2 + \/(0(r ft ))r2 (4.52) 
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, ( , [0 (r h ) /2] 2 [1 + w (r h )] rj-[l- w 2 (r h )] w (r h ) 

w (r h ) = k ^ (4.53) 

r h -[l- w 2 (r h )f /r h - 2 [0 (r h ) /2] 2 [1 + w (r h )f r h - 2V (0 (r h )) r\ 

l(r ) = [l+w(r h )] [0(r fe ) /2\ + V^(nM u 54) 

V r h -[l- w 2 (r h )} 2 /r h - 2 [0 (r h ) /2] 2 [1 + w (r h )] 2 r h - 2V (0 (r h )) r\ 

5' (r h ) = - (2 [w' (r h )} 2 + [r h rf (r h )} 2 ) /r h (4.55) 

on the horizon, and use w (r^) and i] (r/j) as shooting parameters for r^ — 1. The results for 
k = 1 and k = 2 are shown in figs. 4a — d. 

Two solution branches again appear for each k, with the quasi-/c = branch distinguished 
by its weak-gravity v — > limit. The A; = 1 solutions are decribed by < v < 0.356, (1 — 
0.292w 2 ) > w{r h ) > 0.870, -(1 - 0.217*;)*; < 77 (r h ) < -0.304, (0.5 + 1.820*;) < M < 
1.0181, 1.19*; 2 < 5 < 0.3310 (quasi-A; = 0) and < v < 0.331, 0.632 < w (r h ) < 
0.801, -0.826*; > 77 (r h ) > -0.275, 0.9372 < M < 1.0043, 0.5485 > 5 Q > 0.4300. The 
specifications of the k = 2 solutions are < *> < 0.0486, —0.632 < w (77J < —0.518, 
-(1 - 1.62*;)*; > 77 (r h ) > -0.0429, 0.9372 < M < 1.0073, 0.5485 < 5 a < 0.5799 (quasi- 
k = 1) and < v < 0.0475, -0.345 > w(r h ) > -0.489, -0.901*; > r]{r h ) > -0.0417, 
0.9938 < M < 1.0066, 0.5932 > 5 Q > 0.5844. 

We see that these black hole solutions have nontrivial gauge and Higgs field structure 
outside of the horizon. This could not happen if the gauge group were abelian - it relies, as 
discussed in section II, on the non-abelian nature of the field theory being studied here. 



V. Conclusions 

In this paper we have studied spherically symmetric classical solutions to SU(2) non- 
abelian Proca theory and spontaneously broken gauge theory. Our main intent has been 
to exploit a gap in the known no-hair results - they do not necessarily apply to nonlinear 
field theories - to find black hole solutions, which have nontrivial structure outside of the 
horizon, in the most familiar and relevant kinds of field theories . In particular, we have 
presented strong numerical evidence indicating that such field theories, when coupled to 
Einstein gravity, do admit spherically symmetric black hole solutions in which the gauge 
and Higgs fields decay to their asymptotic values exponentially far from the hole, in contrast 
to previous expectations. 
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An important physical question is whether these solutions are stable (although we reem- 
phasize that the no-hair proofs have to do with existence of classical solutions, not with 
their stability). We have not exhaustively studied this question, but it does seem somewhat 
unlikely that they are stable. Our EYMH solutions have a fairly direct relation to sphaleron 
solutions (even though many of them involve strong gravity) and hence are probably un- 
stable. In fact, the arguments of [27-21] and [36] as applied to the smooth and black hole 
solutions of EYM solutions [19,18] seem likely ensure that our solutions are unstable. At 
this time we have only definitively studied this question for the smooth lowest node ENAP 
solutions using linear stability analysis along the lines of [27] which does in fact show these 
solutions to be unstable. It would be of interest to carry out this analysis fully, especially 
in light of the natural suggestion [21] that a no-hair theorem might hold if one demands 
stability. It is also important to develop methods to deal with the latter possibility in a 
general setting. 

Another interesting question is to try to understand at a more fundamental level why 
these smooth and black hole solutions exist. Recently, the authors of [37] have given some 
very interesting plausibility arguments to argue for the existence of the solutions of [19] and 
[18]. These arguments rely on properties of the Yang-Mills configuration space associated 
with the existence of large gauge transformations. However, one can evade the no-hair 
arguments in as simple a theory as two scalar fields so long as the quartic potential is chosen 
judiciously. It would be interesting to see if smooth and black hole solutions can be found 
in such models as this will help to clarify the essential physics. We are presently studying 
this question and will report on it elsewhere. 
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FIGURE CAPTIONS 



1) Two- node regular solutions to Einstein- No n-Abelian-Proca theory. The connection 
(1 + w (r)) and total mass are plotted as functions of radius for a range of vector field 
mass values < fi < /i ma x- Increasing fi in figures (la) and (lb) (and in figure (2) ) 
corresponds to decreasing the value of the radius at which the solution exponentially 
decays to its vacuum value. The k = 2 and quasi-/c = 1 solutions, so named because 
as fi — > they resemble the k — 1 solution of Bartnik and McKinnon [19], bifurcate at 

A 4 = A^max- 

2) One- and two-node black hole solutions to Einstein-Non-Abelian-Proca theory for hori- 
zon radius r\ h = 1. Like the regular case, both odd and even- node solutions ex- 
hibit bifurcation at a maximum value of //. The quasi-A; = solutions resemble the 
Schwarzschild solution w = 1 as fj, — > 0; this is the only class of ENAP solutions with 
a weak-gravity limit. 

3) One- and two-node regular solutions to Einstein- Yang-Mills-Higgs theory for quartic 
coupling A = 1/8 and a range of vacuum expectation values v. Here and in figure (4), 
increasing v corresponds to decreasing the value of the radius at which the solution 
exponentially decays to its vacuum value. As v decreases from its maximal value, the 
quasi-/c = solutions approach a flat space limit which is the well known Yang-Mills 
sphaleron. The other solutions do not admit a weak gravity limit. 

4) One- and two-node black hole solutions to Einstein- Yang-Mills-Higgs theory for quartic 
coupling A = 1/8, a range of vacuum expectation values v and — 1. The quasi- 
k = solution approaches a Schwarzschild black hole metric as v decreases and might 
be heuristically described as a black hole with sphaleron hair. The other solutions do 
not have a weak gravity limit (a limit arbitrarily close to Schwarzschild). 
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